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We present the results of our lattice QCD study of the <?d*Dtt coupling, relevant to the D* — > Dn 
decay. Our computation is made on the gauge field configurations that include Nf = 2 dynamical 
light quarks by using the twisted mass QCD action. From the results obtained at four different lattice 
spacings we were able to take for the first time the continuum limit of this quantity computed on 
the lattice. Our final value, g c = 0.53(3)(3), leads to r(D*+ D°ty + ) = (50 ± 5 ± 6) keV, and is 
in good agreement with the experimental results for the width of the charged D*-meson. 
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I. INTRODUCTION 

The result of the first experimental measure- 
ment of the width of the charged vector D-meson, 
T(D*+) = 96 ± 22 keV |, was a surprise be- 
cause many theoretical predictions of the relevant 
hadronic coupling, gr>*D-m suggested much smaller 
(by about a factor of 3 -r- 4) value for T(D* + ). An 
important novelty this year is that the value of 
T(D* + ) measured at CLEO, has been confirmed 
and made much more accurate by the BaBar col- 
laboration, T{D*+) = (83.5 ± 1.7 ± 1.2) keV @. 

Prior to the CLEO result, the QCD based theo- 
retical estimates of go* Dtt, were made by means of 
various QCD sum rules (QSR). The results of dif- 
ferent QSR techniques were converging to a small 
value of gD'D-K [3n5|, the size of which was chal- 
lenged in ref. [f| where it was argued that it should 
be by about a factor of two larger in order to avoid a 
paradox of non-saturation of the Adler-Weisbcrgcr 
sum rule, which is reasonably well saturated by the 
first couple of states in the case of baryons and light 
mesons. As far as the lattice QCD computations of 
this coupling are concerned, there have been quite 
a few results reported in the static heavy quark 
limit (m c — > oo) @, HI, while only a few calcula- 
tions have been reported for the case of propagat- 
ing charm quark 0, |T(| EH- 1 The first such a 



computation was performed in quenched approxi- 
mation in ref. Q, where the strategy to compute 
9d*Dti has been discussed in detail. It is based on 
the use of the LSZ reduction formula and the ax- 
ial Ward identity, so that gD*D-n can be extracted 
from the computation of the matrix element of the 
light quark axial current, sandwiched by the D* and 
D states. The resulting go-D-m obtained from the 
quenched lattice QCD study (Nf — 0), was large as 
anticipated in Q and compatible with the CLEO 
measurement of F(D* + ). That study has then been 
extended to the unquenched case, i.e. by includ- 
ing the effect of Nf = 2 dynamical quark flavors 
in the QCD vacuum fluctuations Despite the 
fact that the O(a)-improved Wilson quark action 
was used, the results were obtained for relatively 
large pion masses and with a rather poor statis- 
tics at two lattice spacings so that the continuum 
extrapolation was not feasible. In this paper we 
provide the first result of the unquenched compu- 
tation of qp* £>7r by using twisted mass QCD on the 
lattice (mH^i with Nf = 2 dynamical light quarks, 
and at four different lattice spacings. Moreover, our 
computation is made by using considerably lower 
pion masses than in ref. [To| . Thanks to the high 
statistics of the ensembles of gauge field configura- 
tions that we had at our disposal we were able to 
make the continuum extrapolation and obtain the 
physically relevant value of gn* d-k in the continuum 
limit. As we shall see, with respect to the previous 
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1 While D* — > Dtt decay is allowed, B* — > Bit decay is 



kincmatically forbidden. 
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lattice estimates, our gD-D-n coupling is somewhat 
smaller, but still much larger than the predictions 
based on using the QSR techniques. While this 
paper was in writing, an estimate of gD*Dir at a 
single value of the lattice spacing with Nf = 2 + 1 
dynamical O(a)-improved Wilson quarks was made 
in ref. [lH ] which agrees very well with our results 
in the continuum. 

The remainder of this paper is organized as fol- 
lows: we first, in Sec. |TT] define the quantities to 
compute and specify the correlation functions that 
are needed for their extraction; we then in Sec. IIIII 
illustrate the quality of the signals obtained from 
the correlation functions on our lattices and present 
our results obtained for each of the lattice set-ups; 
in Sec. IIVI we discuss the chiral extrapolations and 
compare our result with the other results existing 
in the literature; finally we conclude in Sec. El 



II. DEFINITIONS 

The coupling go-Dir describes the emission of the 
soft P-wave pion off the D* meson, i.e. it is defined 
via the following matrix element 



(D(k)ir(q)\D*(p,\)) = (e x -q) go- 



(1) 



where q = p — k is the pion momentum and A labels 
the polarization state of the vector meson. The 
decay rate of this process is given by 

r(J>* -> Dn) = ^-g-a- g 2 D . D JK\ 3 , (2) 

where C = 1 if the outgoing pion is charged, and 
C = 1/2 if it is neutral, and 



X 1 ^ 2 (m D *,m D ,m 7r ) 



(3) 



where A(a, b, c) = [a 2 — (b — c) 2 ] [a 2 — (6 + c) 2 ], which 
numerically gives \k w -\ = 39.4 MeV and \k v o\ = 
38.3 MeV, i.e. the pion remains very soft. Due 
to T-symmetry and Lorentz invariance, the same 
expression applies to the pion absorption, i.e. to 
Dir —> D* process. 

To compute on* dtt we follow the procedure pro- 
posed in refs. Q, Q, and compute the matrix ele- 
ment of the axial current, = uj^jsd, 



2N e A ■ q u 



{D(k)\A»\D*(p,X)) =2m D *A (q 



A 2 (q 2 )- 



£A • q 



mo + m D , 



?7ln» — m 



D_ n 



(4) 



where ^4o,i,2(<? 2 ) are three independent Lorentz in- 
variant form factors. For q 2 close to ro 2 the reduc- 
tion formula leads to 



(D(k)n(q)\D*(jp,X)) = (D(k)\d^\D*(p,X)) 



which for q 2 = gives, 



gD'Dir 



~17 



4>(0). 



(5) 



(6) 



Since the form factor A (q 2 ) is dominated by pion 
( J p = 0~) in the i-channel, it strongly varies with 
q 2 , as well as with the light quark mass when con- 
sidering the problem for the practical lattice QCD 
calculation. However, since no massless state can 
couple to the axial current, one has 

MO) - ^^MO) mD ;- mD M0),(7) 



2itid 

so that eq. ^ becomes 



go'D-K 



m D 



MO) 



2niD* 



m D , - m D A 2 (0) 



m D Ai(0) 
(8) 



The last expression is suitable for numerical compu- 
tations on the lattice. It reduces to the calculation 
of the following correlation functions 



CW(g;t) = ^(^(0, 0)A V (x, t)P*(y,t s )e 



■iq(x-y)\ 



,(9) 



where the interpolating source operators and 
P$ are placed far from each other at t = and 
t = ts, and create the D*- and D-meson states, 
respectively, at t such that -C t <C ts. The sim- 
plest choice are the local interpolating operators, 
= cy^u and P 5 = ay^d. Note also that q should 
be chosen in such a way that the form factors are 
computed at q 2 = 0, as needed in eq. ©. If the 
vector meson is kept at rest, that means that the 
desired three-momentum should be tuned to 



\q\ 



' D 



2itid* 



(10) 



which is too small a value to be accommodated on 
the periodic lattices explored in the current numeri- 
cal simulations. To remedy this difficulty one of the 
quark propagators, S q (x, 0; U) = (q(x)q(O)), can be 
computed in the rcphascd gauge field configuration, 



U ll (x)^U e Jx) = e ie ^ L U^x) , 



(11) 
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where 6^ — (0,0), and L is the size of the spatial 
side of the cubic box. The propagator 

S S q {x,Q;U) = e^ s ' L S q (x,0;U e ) , (12) 

is equivalent to imposing the twisted boundary con- 
ditions on one of the valence quarks in the consid- 
ered correlation function [l4|, [HI (see also ref . [l6| ) . 
Importantly, however, 9$ in 9 = (9q,9q,9o) can be 
tuned to any value, and therefore also to 



\q\ 



L m 2 D , 



L 



V3 V3 



(13) 



which ensures the extraction of the form factors at 
q 2 = 0. For such a tuned 6q, and by using the 
local interpolating field operators, the correlation 
function ^ reads, 

C^(q;t) = 

(^2 Tr [t/ACO, y)l5S e d (y, x; U)^ u ^S u {x, 0; U) 



x,y 



(14) 



which is very close to what we actually computed 
on the lattice. 2 To relate the above correlation 
function to the relevant form factors we note that 



3 3 



i+3 
(15) 



at large time separations among operators, behaves 
as 



i 

2m 



-™^ 4 x (m D +m D .)^i(0) 



x * D c -En(t s -t) 



2E 



D 



(16) 



where (0\P 5 \D) = Z D , and (0\V t \D*) = efZ D ,. In 
this way we get the first term in eq. ([5]). To reach 
the corrective term in ([5]) we proceed along the lines 
explained in ref. [l(J and compute 



1 3 

C M (q;t)=--J2C t0 (q;t) 



i=l 
3 



mi)< — E 



D 



Cij(q;t) 



(17) 



so that at large separations among the operators, 
one gets 



2 qi m D , A 2 (0) 



C i0 (q;t) _ 
C t3 {q;t) (m D . +m D ) 2 A x (0) ' 



(18) 



where = Qq/L. Hadron masses, m D w, and cou- 
plings, Z D (,), are extracted from the study of the 
large time behavior of the two-point correlation 
functions, namely 



(JXsif; t)V?(0;0)) 



i>0. 



2 cosh[m g * (T/2 - t)] c - m „, T /2 



(^(^^(OjO)) 



t > 0, 



2 cosh[m r) (r/2 - f)] c _ mDT / 2 



(19) 



where T stands for the size of the temporal exten- 
sion of the lattice. Note that in eq. (fT9| we used the 
symmetry of the correlation functions with respect 
to t <-> T - t. 



III. LATTICES USED IN THIS WORK 

Results presented in this work are obtained 
by using the gauge field configurations gener- 
ated by the European Twisted Mass Collaboration 
(ETMC) [13, by simulating the twisted mass QCD 
on the lattice at the maximal twist [l2j]. Nf = 2 
dynamical light quarks have been included in the 
simulations. In tab. U we collect the main informa- 
tion concerning the lattices used in this work. At 
each lattice spacing, simulations have been made 
with several values of the sea quark mass ^ S ea, 
covering the range of the associated pion masses 
m n 6 (280,500) McV [Hj]. Concerning the charm 
quark mass, its value has been discussed in great 
detail in ref. [l8j and it is fixed by requiring the 
agreement between the charmed hadron masses (D, 
D s and r\ c mesons) computed on the lattice with 
their experimentally established values. 

Our first task was to compute the two-point func- 
tions (fT9"|) and extract the accurate values of the 
masses of D q and D* mesons which are necessary 
for determination of 0q in eq. (fTB"]) , 3 In the com- 
putation of the two-point functions ([IT)]) we used 



2 As we shall see in Sec. IIIII the difference with respect to 
eq. (1 1 41 is that in practice we use the extended (smeared) 

source operators. 3 Index "q" in Dj* ' labels the light valence quark, which in 
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source operators composed of local quark fields ip or 
smeared ones tp n . We choose the gaussian smear- 
ing defined via: 



1 + kH 
1 + 6k 



(20) 



where H is the smearing operator defined as [l| 



H * = E (^-w + E>r-i/A-Mj) - (2i) 



with {/"" being the n a times APE smeared link J2(j , 
defined in terms of (n a — 1) times smeared link 
^i,/** anQl i ts surrounding staples V^^° _1 \ 



Proj 



SU(3) 



(1 - a)U (na ~ 1] + -V (n "~ 1] 



For the smearing parameters we choose: 

k = 4, n„ = 30, a = 0.5, n a = 20 , 



(22) 



(23) 



which we kept fixed for all of our lattices. In our 
calculation the quark propagators are computed by 
using stochastic sources, and the correlation func- 
tion are obtained taking advantage of the so called 
one-end trick [l7j | . Furthermore, to improve the sta- 
tistical quality of our results, we use up to eight dif- 
ferent spatial sources. The full list of results is pre- 
sented in tab. UU where we give our values for mn 
and m d* that are totally consistent with our previ- 
ous findings presented in ref. [23| and obtained by 
using the single source propagators. Once we have 
the meson masses we fix the value of 9q according to 
eq. (|13p . and then compute the light quark propa- 
gator with twisted boundary conditions (|12p , which 
is needed to evaluate the three-point functions (fT4|) . 

Instead of the coupling go* Dm it is often conve- 
nient to define g c which is related to gD*D-n via 



2^/m D m Dt 



9c 



(24) 



with the convention that corresponds to / w = 
130.41(3)(20)MeV 0. On the basis of heavy 
quark expansion the coupling g c is expected to scale 
as a constant up to corrections proportional to the 
powers of the inverse heavy quark (meson) mass. 
In our case, we work with the propagating heavy 
charm quark, which is why we added an index "c" 



to avoid ambiguities when comparing g c computed 
here with the values obtained in the static limit Q . 
After combining eq. © with (f2"4"]l we have 



m D * |ffl D 



Mo) 



2^m D m Dt 
G(0) [1 + G'(0)] 



m D , - m D A 2 {0) 



(25) 



where G^{0) arc introduced for shortness. The val- 
ues of the momenta tuned in the twisted boundary 
conditions (fT3"|) are very small and by computing 
G'(0) by using eq. (jTHJ its value is consistent with 
zero for all our lattices. The leading term, G(0), 
instead is extracted accurately from the correlation 
functions (fT5"|) . In fig. [1] wc illustrate the plateaus 
obtained at all of our four lattice spacings and for 
one value of the sea quark mass. Illustrated is the 
ratio 



Ri(t) 



2E D 2m D .Cij(q;t) rt+ED{ts _ t) 
t s >*»o 



}(m D ,+mn.)i 1 (0). (26) 



Corresponding values of g c are given in tab. |TT] 
for all of our lattice set-ups. Fits are made 
on the plateaus corresponding to tg G [8,14] 3 . 9 , 
[11, 14]4.05, [15, 18]4.2, in an obvious notation. 

From the simulations made at /3 = 3.9 with fi q = 
0.004 we were also able to verify the smallness of 
the finite volume effects by comparing the g c values 
obtained on the 24 3 x 48 and on the 32 3 x 64 lattices 
which are completely consistent with each other. 4 

Finally we should also mention that the renor- 
malization factors of the axial current, Za{.9q), 
have been properly accounted for in our compu- 
tation. Their values have been computed non- 
perturbatively in refs. [H|,[2l| and, for convenience, 
they are also given in tab. Q] of the present paper. 



IV. PHYSICAL RESULTS AND 
PHENOMENOLOGY 

A. Getting to the physical value for g c 

The values of the coupling g c , obtained with 
unphysical light quark masses from our lattices, 
should now be related to the physically relevant 



our computations is always kept degenerate in mass with 
the sea quark. 



4 We reiterate that in our computations mass of the sea and 
valence light quarks arc always kept degenerate in mass. 
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FIG. 1: Ratio Ri(t), defined in eg. i(26|) [see also 
eqs. i)J5U6D /. as a function of time, with two sources 
of Dg and D q mesons fixed. Illustration is provided at 
each of the four lattice spacings considered and with the 
light quark mass u, q = 0.008, 0.004, 0.008, and 0.002, 
for f3 = 3.80, 3.90, 4.05, and 4.20, respectively. Also in- 
dicated is the plateau region at which a fit to a constant 
has been made, Ri(t) — > (nriD q + rriD*)Ai(0). 



coupling g c . To do so we will make the contin- 
uum and the chiral extrapolations simultaneously 
by employing three different chiral extrapolation 
formulas. Concerning the continuum extrapolation 
we will rely on the simple linear extrapolation in 
a 2 , which turns out to be adequate for describ- 
ing our data. Concerning the chiral extrapolation 
we will use the simple linear extrapolation in the 
light quark mass (or, equivalently, in the pion mass 
squared), and the two formulas that include the chi- 
ral logarithmic corrections. In other words we will 
use: 



Linear extrapolation: 



g c =«o 1 + K + &— , (27) 



where on the left hand side are the values we 
obtain on our lattices for a given light quark 
mass (pion mass), in which the valence and 
sea quark mass are the same. On the right 
hand side, <?o stands for the same coupling 
obtained in the chiral limit, and a is the slope 
parameter. Notice also that we divide a 2 by 
the value of the lattice spacing obtained at 
f3 = 3.9, so that the /3 parameter is a measure 
of the deviation of our continuum result form 
the values obtained at a = 0.085(3) fm. 

Chiral Perturbation Theory (ChPT-I): 



.9c = .9o 



4g 2 

(4tt/) 2 



loj 



; m 2 + am- 



'3.9 



(28) 



is the formula that has been derived many 
times in the heavy meson chiral perturba- 
tion theory (HMChPT) HH, an effective the- 
ory constructed from the combined heavy 
quark symmetry with the spontaneous chi- 
ral symmetry breaking pattern in order to 
describe the mesons consisting of heavy and 
light quarks [H,[27|. The above formula (gSJ) 
is obtained at the next-to-leading order in chi- 
ral expansion and to the leading order in the 
heavy quark expansion (i.e. in the static limit 
of QCD). Although we are working with the 
charm heavy quark, for which the above for- 
mula is not fully adequate, we will use it as a 
way to guide our extrapolation to the physi- 
cal coupling g c . The constant / in the above 
formula stands for the pion decay constant 
in the chiral limit, which we take to be / = 
120 MeV. In fig. [2] we show the result of ex- 
trapolation guided by eq. (j28j) . For the light 
quark mass we use the definition m^ s (2 GeV) 
in the same way it was defined (renormal- 
ized) in ref. [llj]. As usual, m 2 = 2B n m q , 
and Bp(2 GeV) = 2.71(15) GeV. 

Chiral Perturbation Theory (ChPT-II): 



.90 



2(1 + 2.9q 2 ; 

(47T/) 2 



m 7T l°g' TO 7r 



"3.9 

(29) 



was derived in ref. [271 ] where it was noted 
that the expression (|28|) was obtained by us- 
ing the pion field directly, and not the axial 
current. In principle, our formula (0 is re- 
lated to the pion field directly. However, as 
we mentioned in Sec. Mil our result is com- 
pletely dominated by the axial form factor 
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FIG. 2: Chiral and continuum extrapolation of our lat- 
tice QCD results for g c given in by using eq. 



dence of the heavy-light hadronic quantities, and 
it is even less clear whether or not the above ex- 
pressions could be applied to the case of mesons 
with (not so heavy) charm quark. In other words, 
none of the above formulas is exactly applicable 
to our situation and we will use all of them. The 
spread of results associated with each of the formu- 
las (|27I28|29[) could then be used to estimate the 
systematic uncertainty due to the chiral extrapo- 
lation. That error will leave enough room for the 
future improvement of g c (i.e. of go^On) computed 
on the lattice. 

In more detail, from the fit of our data to the 
above formulas we obtain: 



Linear 



.g = 0.56(3),a 



0.6(2) GcV" 2 , P -- 
0.562(28) . 



-0.12(4), 
(30) 



Ai(0), and therefore the formula (|29|). which 
differs from (|2"5| by the vertex tadpole di- 
agram, should be used in chiral extrapola- 
tion too. Notice that in both cqs. (|28I29|) 
the /^-dependence of the log term and of the 
counter-term coefficient a are not written ex- 
plicitly as they cancel in the sum. 



ChPT - I 

.go = 0.52(2),a = 0.3(2) GeV" 2 , /3 = -0.12(4), 
=> gf iys - = 0.536(25) , (31) 

ChPT - II 
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.90 = 0.46(2),a = -0.16(25) GeV" 2 , 
=> ,9 c phys - = 0.500(24) . 

As our final result we quote 

,9 C = 0.54(3)(^ 2 ), 



0.02 0.03 0.04 0.05 0.06 



m q [GeV] 

FIG. 3: Same as in fig. [2| but using eq. if29)) . 

We stress again that the above expressions are de- 
rived in the static heavy quark limit. It is not clear 
to what extent the chiral logarithmic corrections 
are relevant in describing the light quark depen- 



-0.13(5), 
(32) 

(33) 



where the errors refer to the last digit. The first 
error is statistical and is to be understood as gaus- 
sian, while the second one is uniform and accounts 
for the spread of values obtained by using various 
fit procedures to reach the physical limit. 



B. Phenomenology 



The above coupling (|33[) can now be cast into the 
form in which it is usually discussed in phenomenol- 
ogy by using eq. (|2~4"|) . We obtain, 

go-D, = 9 D , +1 ? V+ = 15.9(7)(± 2 ) , (34) 

which then leads us to our estimate for the decay 
width ©, namely, 



T(D* 



51 ± 5; 



keV. 



(35) 
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The width of the charged D*-meson can be calcu- 
lated as 



T(D* 



=r(zr+ ^^°tt + ) + t{d*+ 
+ r(D*+ -> l>+ 7 ), 



(36) 



as these are the only kinematically available decay 
channels. The last mode, in particular, is a tiny 
fraction of the full decay width and is known ex- 
perimentally, B(D*+ -> D+j) = 0.016(4) 0, so 
that by using our value for g c we finally have 



r(D*+) = (76±7t? ) keV ; 



(37) 



which agrees favorably with the experimental 
values, T(D*+) CLEO = 96 ± 22 keV [H], and 
r(Z?*+) BaBar = (83.5 ± 1.7 ± 1.2) kcV 0. 

Furthermore, we can make a prediction of the 
width of the neutral D* meson, 

T{D*°) = T(D*° -> D-n+) +T(D*° -> DV) 
+T(D*° D° 7 ) , (38) 

if we use the experimental information B(D*° — > 
£>°7) = 0.381(29) [13. Knowing that L>*° -> 
D~ir + is kinematically forbidden, and assuming the 
isospin symmetry to relate our g c to the decay to a 
neutral pion, we obtain 

r(u*)[i - b(d*° ^dS)] = ™ DoKol ) 2 9 2 c 

12tt m D ,o 
= (53 ± 5tr) kcV , (39) 

where \k' o\ — 42.3 McV, and the prime is used 
to distinguish this case, in which all mesons are 
neutral. Our result Q39p is much smaller than 
the relati vely old experimental bound, T(D*°) < 
2.1 MeV [If. 



Comparison with results available in the 
literature 



g c = 0.71(7) was obtained [T(|. However, all these 
results were obtained by working directly with large 
pion masses. Very recently the result of the lattice 
QCD study with C(a)-improvcd Wilson quarks and 
with JVf = 2+1 dynamical light flavors was reported 
in ref. [13j , Although working at the single lattice 
spacing the direct computations were made for a 
broad range of pion masses down to 300 McV. The 
authors fit their result linearly and quadratically to 
extrapolate to the physical pion mass, without rely- 
ing on the expressions derived in HMChPT. Their 
final result, g c = 0.55(6), agrees well with ours, as 
well as with the experimental ones. Current sit- 
uation concerning the estimate of g c is shown in 
fig. 01 where we observe a large difference with re- 
spect to the results obtained by using QSR. Ways 
to remedy this difficulty have been proposed in the 
literature [29| but a real reason for this failure has 
not been fully understood. Here we only focused on 
the computations closely related to QCD. A num- 
ber of results based on quark model calculations 
have been reported in the literature too. For an 
exhaustive list of references see @, |3(J. Extract- 
ing this coupling from the residuum of the semilep- 
tonic form factor was shown to be very difficult be- 
cause the relevant D — > n, B — > n, or D — > K 
semileptonic form factor is a strongly varying func- 
tions at large q 2 -region and the extraction of g c (or 
gb) significantly depends on the parameterization 
used to describe the g 2 -dcpendence of the form fac- 
tor MM- 



V. SUMMARY AND PERSPECTIVES 

In this paper we report on the first lattice QCD 
result for the go* Di? coupling obtained in the con- 
tinuum limit. By using twisted mass QCD on 
the lattice at four different lattice spacings and by 
working with Nf = 2 dynamical light quark flavors, 
we were able to compute this coupling for several 
values of the pion mass, covering the range between 
280 McV and 500 MeV. Our final value is 



To compare our result (|33p with the existing ex- 
perimental ones we use the expression (f36|) and 
from the width reported by CLEO [l| we obtain 
g c = 0.60(6). From the new BaBar result @], in- 
stead, we get g c = 0.56(1). The lattice QCD results 
reported so far were obtained by using the 0(a)- 
improvcd Wilson quarks at single lattice spacing. 
The result obtained in the quenched approxima- 
tion, g c = 0.67(8) (ig) 0, was later confirmed in 
the unqucnchcd study with Nf = 2 dynamical light 
flavors, at nearly the same value of the lattice spac- 
ing, g c — 0.66(11) [lo| . At a slightly finer lattice, 



g c = 0.53(3) (3), & g n , +lf 



15.8(7)(3), 



where we symmetrized the second error bar. That 
value leads to the theoretical prediction of 



T{D* 



(50 + 5 + 6) keV. (40) 



Our value is consistent with the result for g c (or 
g D »+jjO +) extracted from the experimentally mea- 
sured width of the charged vector meson, T(D* ± ), 
after using the isospin symmetry to relate the de- 
cays to a charged and to a neutral pion in the fi- 
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nal state. Furthermore, we were able to predict 
T(D*°) = (53±5±7) keV. 

The second error in all of the above results re- 
flects our estimate of the uncertainty arising from 
the chiral extrapolation, i.e. the way one can relate 
the coupling g c directly computed on the lattice to 
the one corresponding to the physical pion mass. 
That error can be reduced in a new generation of 
lattice QCD simulations performed closer to the re- 
alistic physical situation, namely the lighter pions 
and by including the effect of TVf = 2 + 1 + 1 dynam- 
ical quarks in the QCD gauge field configurations. 
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p 


3.8 


3.9 


3.9 


4.05 


4.2 


4.2 


L 3 x T 


24 3 x 48 


24 3 x 48 


32 3 x 64 


32 :! x 64 


32 3 x 64 


48 3 x 96 


# meas. 


240 


240 


144 


144 


144 


96 


/-^scal 


0.0080 


0.0040 


0.0030 


0.0030 


0.0065 


0.0020 


A l sca2 


0.0110 


0.0064 


0.0040 


0.0060 






/-^sca3 




0.0085 




0.0080 






/-^sca4 




0.0100 










a [fm] 


0.098(3) 


0.085(3) 


0.085(3) 


0.067(2) 


0.054(1) 


0.054(1) 


Za( 5o 2 ) [21, 21] 


0.746(11) 


0.746(6) 


0.746(6) 


0.772(6) 


0.780(6) 


0.780(6) 




0.2331(82) 


0.2150(75) 


0.2150(75) 


0.1849(65) 


0.1566(55) 


0.1566(55) 



TABLE I: Lattice ensembles used in this work with the indicated number of gauge field configurations. Lattice spacing is 
fixed by using the Sommer parameter ro/a |22| , with tq = 0.440(12) fm fixed by matching /„■ obtained on the lattice with 
its physical value (c.f. rcf. [IS}). Quark masses are given in lattice units. 
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ft 
p 


/ 




MS/o PpV^ 

( /tg \jKi V ^ 


m D q 




9c 


3.80 


24 


/*seal 


0.0398(11) 


1.73(5) 


1.95(6) 


0.536(36) 






/isca2 


0.0547(15) 


1.75(5) 


1.98(6) 


0.524(26) 


3.90 


24 




0.0216(5) 


1.75(5) 


1.96(5) 


0.559(33) 






/isea2 


0.0345(8) 


1.77(5) 


1.98(5) 


0.564(24) 






/^sea3 


0.0458(11) 


1.78(5) 


1.99(5) 


0.589(18) 






/^sca4 


0.0539(13) 


1.78(5) 


1.99(5) 


0.588(13) 


3.90 


32 


/^scal 


0.0162(4) 


1.74(5) 


1.94(5) 


0.509(24) 






/^sea2 


0.0216(5) 


1.75(5) 


1.95(5) 


0.543(37) 


4.05 


32 


/-^scal 


0.0249(7) 


1.81(4) 


2.01(4) 


0.528(29) 






Msea2 


0.0374(10) 


1.83(4) 


2.03(5) 


0.585(22) 






/^sca3 


0.0499(13) 


1.84(4) 


2.04(5) 


0.575(22) 


4.20 


32 


/^scal 


0.049(2) 


1.91(4) 


2.11(5) 


0.620(18) 


4.20 


48 


/^scal 


0.0150(7) 


1.87(4) 


2.06(4) 


0.540(35) 



TABLE II: D q and D* meson masses obtained with the charm quark mass Bxed to its physical value and the light 
quark mass equal to the mass of the sea quark, specified in tab. [I] and for various light quark masses. The values 
for g c are obtained as described in the text. All masses are given in physical units [GeV], and g c is dimensionless. 



